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Abstract 

Previous studies have shown that extrasolar Earth-like planets in close-in habit- 
able zones around M-stars are weakly protected against galactic cosmic rays (GCRs), 
leading to a strongly increased particle flux to the top of the planetary atmosphere. 
Two main effects were held responsible for the weak shielding of such an exoplanet: 
(a) For a close-in planet, the planetary magnetic moment is strongly reduced by 
tidal locking. Therefore, such a close-in extrasolar planet is not protected by an 
extended magnetosphere. (b) The small orbital distance of the planet exposes it to 
a much denser stellar wind than that prevailing at larger orbital distances. This 
dense stellar wind leads to additional compression of the magnetosphere, which can 
further reduce the shielding efficiency against GCRs. In this work, we analyse and 
compare the effect of (a) and (b), showing that the stellar wind variation with orbital 
distance has little influence on the cosmic ray shielding. Instead, the weak shielding 
of M star planets can be attributed to their small magnetic moment. We further 
analyse how the planetary mass and composition influence the planetary magnetic 
moment, and thus modify the cosmic ray shielding efficiency. We show that more 
massive planets are not necessarily better protected against galactic cosmic rays, 
but that the planetary bulk composition can play an important role. 
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1 Introduction 



One of the many fascinating questions in the field of exoplanet studies is 
the search for habitable worlds. Because of their relatively small mass, low 
luminosity, long lifetime and large abundance in the galaxy, M dwarfs are 
sometim es suggested as prime targets in se arches for terrestrial habitable 



planets (iTarter et al.l . 120071 : IScalo et all 120071 ). The detection of planets with 
M < IOMe, i.e. with a mass smaller than 10 terrestrial masses (t he upper 



limit for "super- Earths" in the classification of [Valencia et al.l . 120071 ) recently 



became possible with current advances in instrumentation and analysis. The 
first super-Earth detecte d was GJ 876d, a planet with ~7.5Me orbiting an M 



star (IRivera et al.l . |2005| ). Thereafter, other low- mass planets have been dis- 



covered: OGLE-2005-BLG- 390Lb, a planet with 5.5Mk in an orbit of 5 AU 



was found by microlensing (IBeaulieu et al.l . 120061 ). and HD 69830b, a planet 



with a minimum mass of 10.2Mf. in an orbi t of 0.08 AU was detected by the 
radial velocity technique (iLovis et al.l . l2006l ). More recently, two super-Earth 
planets were discovered around the M3 dwarf Gl 581, with mini mum masses 



of 5.0 Me and 7.7 Me, and se mimajor axes of 0.07 and 0.25 AU (lUdry et al. 



20071 ). and lRibas et al.l (120081 ) report the possible detection of a 5.0 Me planet 
around the M-type star GJ 436. In total, 6 planets with masses < IOMe are 
known today, 5 of which are located in an orbit around an M stan ^ I. The least 
massive planet has a mass of 5.0 Me. While the number of super-Earth plan- 
ets around M stars is currently still limited, the detection of cold debris disks 
around M stars seems to indicate that planets may be a s frequent for M stars 
as they are for F, G and K stars (ILestrade et al.l . l2006l ). 



For M stars the habitable zone (defined by the range of orbital distances over 



which liquid water is possible on the planetary surface, see e.g. iKasting et al. 



19931 ) is much closer to the star (depending on stellar mass, but typically 
<0.3 AU). Such close-in distances pose additional problems and constraints 
to habitability. For example, there is the possibility that the atmosphere might 
collapse on the night side of the plane t, although it seems that this p roblem is 



less severe than previously assumed (jJoshi et al.l . 119971 : iJoshil . |2003| ) . Intense 



stellar flares are common, especially for young M stars, but the radiation can 
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be shielded b y the planetary atmo sphere (IHeath et al.l . Il999l ). Both direct 
UV radiation (IBuccino et al.l. 120071) an d indirect UV radiation generated by 



energetic particles (jSmith et al.l . |2004| ) can have important consequences for 



biogenic processes on M-star planets. Another potential problem for M star 
planet habitability is that one expects a large number of coronal mass ejections 
(CMEs) on active M stars. The related interaction of the dense CME plasma 
flux with the atmosphere/magnetosphere environment of the exposed planets 
during the active stage of the stellar evolution could be strong enough to erode 
the atmospher e or the planetary water inventory via non-therra al atmospheric 
loss processes (iKhodachenko et al.l . 120071 : iLammer et al.l . l2007l ). 



The efficient loss of water from the planet could also lead to a breakdown of 
plate tectonics. If the majority of the planets' water is lost, the lithosphere will 
not be sufficiently deformable, and subduction of the crust cannot occur. The 
recycling of the crust through plate tectonics keeps the crust thin. If the crust 
is too thick, the lithospheric plate comprising the crust will be too buoyant to 
be subducted. A thin crust and a wet planet seem to be very important for 
plate tectonics to operate. Plate tectonics, however, help to cool the interior 
efficiently which is required to keep an Earth-lik e thermally driven rn agnetic 



dynamo in operation over geological time spans (ILammer et al.l . |2008| ). Even 



in the case of M star planets with a magnetic dynamo, the magnetic moment 
is not likely to be large. Terrestrial planets in the habitable zones of M-stars 
can be expected to be ti dally locked, which leads to a reduction in the plane- 
tary magnetic moment (iGrieBmeier et al.l . l2005bl ). For these reasons, one can 
assume that M-star planets within the habitable zones are only weakly pro- 
tected by their magnetospheres, resulting in a higher cosmic ray exposure. In 
addition, the small orbital distances of these planets expose it to much denser 
stellar winds than at larger orbital distances. The dense stellar winds lead 
to a further compression of the magnetospheres, which can in principle again 
reduce the shielding efficiency against GCRs. 



In this work we extend the studies of lGrieBmeier et al.l (l2005bl ) and lGrieBmeier et al, 
( l2005al ) by separately quantifying and comparing these two effects (effect of 
enhanced stellar wind ram pressure due to the smaller orbital distance vs. in- 
fluence of a reduced planetary magnetic dipole moment). We show that the 
stellar wind variation with orbital distance has little influence on the shielding 
efficienc y against galactic cosmic rays. The weak shielding found for M star 



planets ( GrieBmeier et al. . 2005bl 



a 



thus has to be attributed to their reduced 
magnetic moment. Given that the magnetic moment is the key parameter, we 
further analyse the influence of the planetary mass and its composition on the 
planetary magnetic moment, and thus on the cosmic ray shielding efficiency. 



This paper is organised as follows: The planetary parameters governing the size 
of the magnetosphere are presented in Section [2l We show that planets within 
the habitable zone of M stars are tidally locked, which leads to a reduced 
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planetary magnetic dipole moment (Section I2.2p . The stellar wind parameters 
at close orbital distances are discussed in Section 12. 3[ The influence of tidal 
locking and stellar wind on the size of the planetary magnetosphere is briefly 
mentioned in Section I2.4[ In Section I2.5[ we discuss the dependence of the 
cosmic ray flux outside the planetary magnetosphere on the orbital distance. 
The cosmic ray model {inside the magnetosphere) is explained in Section [31 
In Section |H the resulting flux of galactic cosmic rays to the atmospheres of 
different extrasolar planets is discussed. We study the effect of orbital distance 
via the stellar wind parameters (Section [4. ip . tidal locking (Section [4. 2p . plan- 
etary mass and planetary composition (Section l4.3l) . Potential implications are 
briefly discussed in Section |5l Section [6] closes with some concluding remarks. 



2 The planetary situation 



2.1 Tidal locking 



Tidal locking (i.e. the synchronisation of the planetary rotation period to 
the planetary orbital period) occurs for planets orbiting very close to their 
parent star. The relevance of this process for terrestrial extrasolar planets is 
discussed in this section . This is an update to the corresponding section of 



Griefimeier et al.l fl2nn5bl ). 



A planet with angular velocity c^i at t = (i.e. after formation) will gradually 
lose angular momentum, until the angular velocity reaches a constant value 
LJ{ at t = Tsync- For a planet with a mass of Mp and radius of Rp around a 



star of mass M^, Tsynr can be expressed as (e.g. Murray and Dermottl . Il999 



Griefimeier et al.l . 120071 . Appendix B): 

4 „, / ^ 



sync - Q-^PyQMp 




The parameters a, Q'^, oji and ojf required to calculate the timescale for tidal 
locking of a terrestrial planet are presented further below. 

Eq. ([T]) shows that the timescale for tidal locking strongly depends on the 
distance (rgync oc d^). Thus, a planet in a close-in orbit around its central star 
is subject to strong tidal interaction, leading to gravitational locking on a very 
short timescale. Other important factors for the tidal locking timescale are the 
stellar mass and the planetary structure. 
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2.1.1 Structure parameter a 



The constant a depends on the internal mass distribution within the planet. 
It is defined by a = I/lMpR"^), where I is the planetary moment of inertia. 
For a sphere of homogeneous density, a is equal to 2/5. For planets, generally 
a < 2/5. For the Earth, the structure parameter a is given by a = 1/3 
( iGoldreich and Soterl . Il966l ). In the following, this value will be used. 



2.1.2 Tidal dissipation factor Q' 



is the modified Q- value of the planet. It is defined by (IMurray and Dermott 
1999h 

Q', = (2) 



2h 



2,P 



where /c2,p is the Love number of the planet. Qp is the planetary tidal dissi- 
pation facto r (the large r it is, the smaller is t he tid al dissipation), defined by 
MacDonaldl (Il964j ) and IGoldreich and Soterl (119661 ). 



We u se /c2,p = 0.3 for the Earth (iMacDonaldl . |1964J : iMurray and Dermott 
19991 ). This value will also be used for the "small super-Earth" and the "Ocean 
Planet" case. For the Earth, a value of Qp ~ 12 can be determined from 



the measured secular acceleration of the moon ( Goldreich and Soterl . Il966 



Murray and Dermottl . Il999l ). This value is relatively small when compared 
to other terrestrial planets, where Qp is typically of the order of 10^. This 
is probably due to the fact tha t , for t he Earth, much energy is diss i pated 



in the shallow seas (iMacDonaldl . 1 19641 : iHubbardl . |1984J : iKasting et al 



1993 



Murray and Dermott . 19991 ). Inth e past, when t he continents were ioined, th e 



value of Qp was probably larger (jPealel . Il999l : iMurray and Dermottl . Il999h . 
For the case of the "small super-Earth" and the "Ocean Planet" (see below) 
Qp is set to the value usua lly assumed for Venus and Mercury {Qp ^ 100 



Murray and Dermottl . Il999l ). 



From fc2,p and Qp, the required value of Q'^ can be obtained using eq. ([2]). We 
find Qp = 60 for Earth and Qp = 500 for both the "small super-Earth" and 
"Ocean Planet" cases. 



2.1.3 Initial rotation rate uji 

The initi al rotation rate uj\ of a ter restrial planet is a poorly known quantity 



see, e.g. ICorreia and Laskarl . |2003| ). It certainly depends on the details of the 



planetary formation and can be strongly infiuenced by processes like migration 
or impacts. Therefore, two limits for uJi are considered: 



A relatively high initial rotation rate as suggested for the early Earth- Moon 
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system, with oji = 1.83 cje (with the current rotation period of the Earth 
ujk = 7.27-10"^ s~^) corresponding to a length of day of 13.1 h (IMacDonald . 



19641 ) 



A lower rotation rate with cui = 0.80 cue corresponding to a day of 30 h. 



Note that a primordial rotation period of the order of 10 h is consistent with 
the relation between the planetary angular n i omeii tum density and planetary 
mass observed in the solar system (IHubbardl . Il984j . Chapter 4) . 



2.1.4 Final rotation rate 



As far as eq. ([T]) is concerned, the final rota tion rate Uf can b e neglected, at 
least for the planets of interest in this work (IGriefimeierl . |2006| ) . 



2.1.5 Tidal locking: Results for exoplanets 

Eq. ([1]) allows us to calculate the timescale for tidal locking for terrestrial exo- 
planets. Subsequently, the planets can be classified as either "tidally locked" , 
"potentially tidally locked" or "unlocked". The upper and lower boundaries 
of the potentially locked region are determined by the conditions Xsync = 100 
Myr and Tgy^c = 10 Gyr, respectively. To account for the uncertainty of the 
initial rotation rate cu;, the lower boundary is calculated with cu; = 1.83 cje 
(i.e. a rotation period of 13.1 h), and the upper boundary is calculated with 
Ui = 0.80 We (corresponding to a rotation period of 30 h). This increases the 
of the "potentially locked" region. 



In Figure [T] (updated from iGrieBmeier et al.l . l2005bl ). the grey-shaded area 



gives the loca t ion o f the continuously habitable zone (GHZ) according to 



Kasting et al.l (119931 ). The lines delimit different parameter regimes for an 
Earth-like planet: All planets to the right of the dashed line are considered 
as tidally locked while planets to the left of the solid line are supposed to be 
freely rotating. One can see that all Earth-like planets inside the GHZ of M 
stars can be considered to be tidally locked. 

In addition to the current Earth, two additional configurations are examined: 

• The "small super-Earth" case (described in more detail in Section Hj) is 
depicted in Fig. [21 It can be seen that the larger value of Q'^ slightly decreases 
the "tidally locked" regime. 

• The "Ocean Planet" case (described in more detail in Section Hj) is depicted 
in Fig. [Sj with similar results than for the "small super-Earth" case. 

As a typical test terrestrial exoplanet at an orbital distance of 0.2 AU 

around a star of 0.5 Mq will be studied. Figs. [T] to [31 clearly show that such 
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a planet, regardless of its precise mass, radius, and composition, is very likely 
to be tidally locked. The reduced rotation rate of tidally locked planets is 
supposed to have important consequences for the planetary magnetic dipole 
moment, which is discussed in Section fI72[ 




Fig. 1. Tidally locked (left) vs. freely rotating (right) regime for Earth-like plan- 
ets as a function of orbital distance d and mass of the host star. The shaded 
area gives the lo c ation o f the continuously habitable zone (CHZ). Updated from 



Griefimeier et all (|2005bl ). 



2.2 Planetary magnetic moment 



The p lanetary magnetic dipole moment can be estimated from different scaling 



laws (IGriefimeier et al.l . l2005bl ). It can be shown that for the slowly rotating 



tidally locked planets, the magnetic moment is much smaller than for freely 
rotating planets. For an Earth-like planet in an orbit of 0.2 AU around a star 
with 0.5 stellar masses, the magnetic moment is expected to lie in the range 
0.02A1e < -M < 0.15A^F,i where M .p, is the value of Earth's current magnetic 
moment (IGrieBmeier et al.l . l2005bl ). In the following, the maximum value of 
O.ISA^E will be adopted to obtain a lower limit for the cosmic ray flux to the 
atmosphere. Magnetic moments resulting for a few other configurations are 
given in Table [H 



7 




Fig. 2. Tidally locked (left) vs. freely rotating (right) regime for a "small su- 
per-Earth" as a function of orbital distance d and mass of the host star. The 
shaded area gives the location of the continuously habitable zone (CHZ). 




J[AU] 



Fig. 3. Tidally locked (left) vs. freely rotating (right) regime for an "Ocean Planet" 
as a function of orbital distance d and mass of the host star. The shaded area 
gives the location of the continuously habitable zone (CHZ). 

2.3 Stellar wind 



It is known that at the close orbital distances of M-star habitable zones the 
stellar wind has not yet reached the quasi-asymptotic velocity regime. Because 
of the low stellar wind velocity, the planetary magnetosphere is less strongly 
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compressed at such distances than one would expect from stellar wind models 
with a constant velocity. To capture this behaviour realistically and to cor- 
rectly describe the flux of galactic cosmic rays into the atmospheres of close-in 
exoplanets, we require a model with a radially dependent stellar wind velocity. 

For slowly rotating stars, the stellar wind may be descr i bed b y the solution of 
the purely hydrodynamic, isothermal model of iParkerl (119581 ). In this model, 
the interplay between stellar gravitation and pressure gradients lea ds to a 
supersonic flow fo r sufficiently l arge substellar distances d (see e.g. iParkerl . 



19581 : iMann et al.l . Il999l : iProlsd . |200J). In accordance with the observations. 



this model describes a solar wind with low velocity and large acceleration near 
the sun, whereas at larger distances the velocity is large and the acceleration 
strongly decreases. 

The procedure to obtain the stellar wind velocity Ves{d, M^, i?*) and density 
n{d, M^, i?^) at the location of an exoplanet (i.e. at distance d) for a host star 
of given mass M^, and radius consists of the following: 

(1) The stellar mass loss rate is assumed to be proportional to its surface 
area: 

M. = MeS. (3) 

where the subscripts and denote solar and stellar properties, respec- 
tively. 

(2) A Parker-like stellar wind model is used to find n{d) and v{d) as a function 
of the distance to the star. The coronal temperature Tcomna is adjusted 
until the stellar wind velocity at 1 AU corresponds to the value which 
was obtained in step 1. With this value of Tcorona, v can be determined for 
any value of d by solving Parker's wind equation. Thus, v{d) is obtained. 
Because of = And"^ n{d)v{d)m (where m is the proton mass), the 
density n{d) is obtained by dividing the stellar mass loss obtained in 
step 1 by 4:7id'^ v{d)m, where v{d) was obtained in step 2. 
The effective stellar wind velocity relative to the planet Ves is calculated 
from V (as obtained in step 2 above) and the planetary orbital velocity 



(3) 



(4) 



This procedure is similar to that of iGrieBmeier et al.l (j2005bl ). with the differ- 
ence that here we keep the stellar age fixed (at 4.6 Gyr) to compare stars of 
identical ages. 



2.4 Magnetospheric model 



The r nagnetosphere is rnodelled as a cyhnder topp ed by a half-sphere (IVoigt 



198ll ; IStadelmannl . l2005l ; IGrieBmeier et al.l . l2005bl ). A closed magnetosphere 
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is assumed, i.e. field lines are not allowed to pass through the magnetopause. 
Within this model, the magnetic field is defined for any point inside the mag- 
netosphere as soon as the planetary magnetic moment and the size of the 
magnetosphere are prescribed. The size of the magnetosphere is characterised 
by the magnetopause standoff distance Rg, i.e. the extent of the magneto- 
sphere along the line connecting the star and the planet. Rg can be obtained 
from the pressure equilibrium at the substellar point. This pressure balance 
includes the stellar wind ram pressure, the stellar wind thermal pressure of 
electrons and protons, and the planetary magnetic field pressure: 

W„ + 2„*,^ = ^||J. (4) 

Here, /o = 1.16 is the form factor of the magnetosphere and inclu des the 



magnetic field caused by the currents fiowing on the magnetopause (IVoigt 



19951 : iGrieBmeier et al.l . |200J). Ai is the planetary magnetic dipole moment 



as discussed in section 

For a given planetary orbital distance d, only the magnetospheric magnetic 
pressure is a function of the distance to the planet, while the other contribu- 
tions are constant. Thus, from the pressure equilibrium eq. (jl]) the standoff 
distance Rg is found to be: 



R.<s 



Svr^ [ranvl^ + 2 nkBT) 



(5) 



Table [T] lists the standoff distances Rg (in planetary radii Rp) for different 
planetary configurations as well as the magnetic field strength at the location 
of the magnetopause, B{Rs). 



2. 5 Cosmic ray flux outside the magnetosphere as a function of orbital dis- 
tance 



In the solar system, the cosmic ray fiux is known to increase with increasing 
heliocentric distance. This is attributed to effects such as diffusion, convec- 
tion, adiabatic deceleration as well a s gradient and curvature drifts (see, e.g. 



Kallenrode. 2000: Heber et al 



20061 ). The theoretical basis for this approach 



is described by iFichtnerl ( 1200 ll . section 2.3). In this section, we study how 
strong such effects are in the inner heliosphere to see how much smaller the 
cosmic ray flux outside the magnetosphere of a close-in exoplanet (e.g. at 0.2 
AU) may be compared to one in an Earth-like orbit. The penetration of these 
particles through the planetary magnetosphere is described in section |H 



The radial gradient of the cosmic ray intensity J depends on the heliocentric 
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Fig. 4. Relative cosmic ray flux in the heliosphere as a function of heliocentric 
distance. 



distance d. Between 1 and 42 AU, iMcDonald et al.l (119921 ) approximate it by: 

(6) 



1 dJ _ ^ fde 
J dd ^ \ d , 



During solar minimum, a ^ 0.7 for cosmic ray protons (the exact value de- 
pends on the particle energy), and Go 
Earth is denoted hj de = I AU. 



0.2 AU ^. The orbital distance of the 



Eq. ([HD is equivalent to a relative cosmic ray flux of: 



Jrel{d) = Joexp 



a 



(7) 



where Jo is determined by the condition Jrei(c?e) = 1- Figure H] shows Jrci{d). 

With this, we can compare the cosmic ray flux at different distances. We find 
that J(0.2AU) ^ 0.8 • J(l.OAU), and J(40AU) ^ 4 • J(l.OAU). For distances 
below 1.0 AU, the particle flux varies only weakly, so that we can neglect this 
effect when calculating fluxes of cosmic rays into the atmospheres of close-in 
planets. 

Note that the radial gradient could be different for other stars than in the solar 
system. For a low mass M star with a lower stellar wind density, we would 
expect the radial gradient of the cosmic ray intensity to be smaller, whereas it 
is likely to be larger for a star with a denser and faster stellar wind (as would be 



the case, e.g. for a younger star, see lWood et al.l . |2002| . l2005l : iGriefimeier et al. 
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2005bl : iHolzwarth and Jardind . 120071 ) . While a quantitative computation of the 



cosmic ray gradients around other stars would be interesting, this is outside 
the scope of this paper. This question will be addressed in the future. For 
the moment, we will assume that the cosmic ray modulation around other 
stars is not much larger than that found in the solar system, and we will 
correspondingly neglected it in the remainder of this work. 



3 Cosmic ray protection 



In order to quantify the protection of extrasolar Earth-like planets against 
galactic cosmic rays, the motion of galactic cosmic protons through planetary 
magnetospheres is investigated numerically. 



3.1 Cosmic my calculation 



In order to determine the impact of GCR protons in the energy range 64 
MeV < < 8 GeV on the planetary atmosphere, particle trajectories in 
the magnetosphere are analysed. Because no solution in closed form exists, 
this is only possible through t he numerical integration of many individual 



trajectories (ISmart et al.l . l2000l ). For each particle energy (64 MeV, 128 MeV, 
256 MeV, 512 MeV, 1024 MeV, 2048 MeV, 4096 MeV and 8192 MeV) and for 
each magnetospheric configuration, over 7 million trajectories are calculated, 
which correspond to protons with different starting positions and starting 
velocity directions. The particles are launched from the surface of a sphere 
centred on the planet. The radius of this sphere satisfies the condition r > Rs, 
so that the particles are launched outside the magnetosphere (except for those 
arriving from the tailward direction). 

Once the particles enter the magnetosphere, their motion is influenced by the 
planetary magnetic field. This magnetic field is calculated from the magne- 
tospheric model of Section 12. 4[ Note that Table [T] gives the maximum values 
for the planetary magnetic moments. For this reason, the results represent the 
lower limit for the cosmic ray flux to the atmosphere. For a smaller magnetic 
moment, a larger cosmic ray flux is possible. The trajectories are calculated 
using the numerical leapfrog method. For each energy, all particles are counted 
which reach the atmosphere (described by a spherical shell one hundred kilo- 
metres above the planetary surface, i.e. = Rp + 100 km). The impact of 
particles on the planetary atmosphere is quantified by the energy spectrum, 
which is defined in Section 13. 2[ More details on the numerical calculation of 



the cosmic rays trajectories can be found in IStadelmannI (120051 ). 
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Note that for the planetary system we imphcitly assume an environment 
similar to that of the solar neighbourhood. Planets in other regions of the 
galaxy (e.g. those encountering a different interstellar medium, or those close 
to the galacti c centr e) may be subject to very different cosmic ray fluxes (e.g. 
Scherer et aLl . l2008h . 



3.2 Cosmic ray energy spectrum 



The cosmic ray energy spectrum is determined in the following way: for a given 
particle energy, the fraction of particles reaching the planetary atmosphere is 
registered. This value is compared to the fraction of particles reaching the 
atmosphere of an identical, but unmagnetised planet. The resulting magneto- 
spheric filter function is multiplied with the co smic ray energy spectrum out- 
side the magnetosphere, which w as taken from ISeo et al.l (119941 ). Fig. [5] shows 
the reference energy spectrum of ISeo et al.l (119941 ) as a dash-dotted line. The 
energy spectrum on top of the Earth's atmosphere is shown as a solid line, 
whereas the dashed line corresponds to the cosmic ray energy spectrum for an 
Earth-like exoplanet in an orbit of 0.2 AU around a 0.5 Mq K/M star. 




proton energy [MeV] 



Fi g. 5. Das h -dotte d line: energy spectrum outside the magnetosphere, according 



to 



Seoetal 



Dashed line: energy spectrum of cosmic ray protons impact- 
ing the atmosphere (100 km above the surface) of an Earth-like planet at 0.2 AU 
around a 0.5 Mq K/M star. Solid line: energy spect r um of cosmic rays impacting 



the atmosphere of the Earth. From Griefimeier et al. ( 2005bl ). 
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4 Comparison of GCR fluxes for eight scenarios 



In the following, the flux of cosmic rays at the top of the atmosphere will be 
compared for eight different planetary conflgurations. The energy spectrum is 
determined for the following cases: 

• Case 1: "Unmagnetised case". This is the (reference) cosmic ray energy 
spec trum outside any planetary magnetosphere (interplanetary case, taken 



from lSeo et al.l . ll994j ). At the same time, this corresponds to the flux arriving 
at the top of the atmosphere of a totally unmagnetised planet (i.e. without 
even induced magnetic flelds). This case is valid for any stellar mass. In the 
following flgures, it is represented by a dash-dotted line. 
Case 2: "(Unlocked) Earth case". The case of the magnetosphere of today's 
Earth (i.e. at 1.0 AU of a sun-like star with = I.OMq). The Earth's 
magnetic fleld is represented by a zonal dipole with a Gaussian coeflicient 
of gio = 31 /iT. In all flgures, this case is represented by a solid line. 
Case 3: "Locked exoplanet case". The case of the magnetosphere of an 
extrasolar planet in an orbit of 0.2 AU around an K/M star of 0.5 solar 
masses. Due to tidal locking, the magnetic moment is reduced to 15 % of 
the Earth's magnetic moment (see Section [2^21) . For this situation, a dashed 
line is used. 

Case 4: (Hypothetical) "unlocked exoplanet case". Identical to the "locked 
exoplanet case" (e.g. located at 0.2 AU around a sun-like star with = 
0.5Mq), with the only difference being the magnetic moment: Here we as- 
sume that the planet is not tidally locked and has a strong magnetic moment 
{Ai = I.OA^e)- Note that this case is not realistic (unless the planet was 
brought to this position only recently), but it is instructive to compare the 
influence of the stellar proximity and the small planetary magnetic moment. 
In the flgures, this case will be shown by fllled triangles. 
Case 5: (Hypothetical) "locked Earth case". Identical to the "(unlocked) 
Earth case" (e.g. at 1.0 AU distance from a sun-like star with = 
I.OM0), but assuming that the planet only has a small magnetic moment 
{Ai = 0.15 Ai^). This case can also be seen as a "weakly magnetised Earth 
case" . Similarly to the "unlocked exoplanet case" , this case will be used to 
disentangle the effect of stellar proximity and of small magnetic moment. 
For this case, circles are used. 

Case 6: "Sr aall super-Earth c ase" (6 Earth masses, 1.63 Earth radii, as 



modelled by lLeger et al.l . l2004l ). Similar to the "locked exoplanet case" (case 
3), but for a large terrestrial exoplanet, assuming tidal locking. As was 
mentioned before, the flrst planets of such size have now been detected. 
This case is represented by a double-dashed line. 

Case 7: "Large super-Earth ca se" (10 Earth masses, 1.86 Earth radii, making 



use of the scaling relations of [Valencia et al.l . l2006l ). Similar to the "locked 



exoplanet case" (case 3) or the "Small super-Earth case" (case 6), but for a 
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larger planetary mass, assuming tidal locking. This planet is denoted by a 
triple-dashed line. Together with the "locked exoplanet case" and the "small 
super-Earth case", this will be used to analyse the dependence of cosmic 
ray protection on the planetary mass. 

Case 8: "Ocean plan et case" (6 Earth masses, 2.0 Earth radii, as modelled 



by lLeger et al.l . |2004| ) . Similar to the "Small super-Earth case" (case 6), but 



for an Ocean Planet (i.e. a planet with a significaiit fract ion of water ice, see 



e.g. iKuchnerl . l2003l : iLeger et al.l . |2004 ISelsis et al.l . 120071 ) . For this situation 



a dotted line is used. 



For each of the above cases, the input parameters (the orbital distance d and 
the stellar mass M^,) and the resulting values (the magnetic moment Ai, the 
standoff distance Rg and the magnetic field strength at the magnetopause 
B{Rs)) are given in Tabled! The stellar mass has two functions: It determines 
the rotation rate of tidally locked planets (case 3 and cases 5-8), and it defines 
the stellar wind flux onto the magnetosphere, thus defining the size of the 
magnetosphere via Eq. (3). For case 1 (no magnetosphere), the stellar mass 
does not matter. As noted before, the magnetic moments (and thus also the 
standoff distances) are assumed to have the maximum allowed value. Thus the 
values obtained in this section represent a lower limit for the flux of cosmic 
ray protons to the planetary atmosphere. Table [1] also includes the magnetic 
field strength at the magnetopause. By definition of the magnetic moment, 
B{Rs) oc /iAR~^. When the contribution of the thermal pressure in eq. (jlj) is 

negligible, this results in B{Rs) oc Jnvl^, so that the magnetic field strength 



at the magnetopause is identical e.g for cases 3 and 4, but different for case 5. 



The cosmic ray flux to the planetary atmosphere depends also on the stellar 
age (which is a crucial par ameter for the ste l lar win d velocity and density). 
This effect was discussed in Griefimeier et al. (l2005bl l 



Influence of orbital distance via the stellar wind parameters 



Previously, the cosmic ray flux to the atmosphere of a tidally locked Earth- 
l ike exoplanet around a K /M type star with = O.5M0 has been calculated 
fiGrieBmeier et al.l . l2005bh . ]t was shown that such a planet (corresponding to 
our case 3) will be subject to very high cosmic ray fluxes when compared to the 
Earth (corresponding to our case 2). For particle energies below 200 MeV, the 
cosmic ray flux to the exoplanet was found to be up to one order of magnitude 
higher than on Earth, and for energies above 2 GeV magnetospheric shielding 
is negligible for the e xoplan et case. This can also be seen in Fig. [51 As stated 
by lGrieBmeier et al.l (l2005bl ). two effects can contribute to this high cosmic ray 
flux: (a) the reduced planetary magnetic dipole moment due to tidal locking 
(Section 12.21) and (b) the enhanced stellar wind ram pressure at small orbital 
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Planet 


d 




M 








[AU] 








[nT] 


Case 1: Unmagnetised 


any 


any 





- 


- 


Case 2: (Unlocked) Earth 


1.0 


1.0 


1.0 


9.91 


73 


Case 3: Locked exoplanet 


0.2 


0.5 


0.15 


4.12 


150 


Case 4: (Hyp.) unlocked exoplanet 


0.2 


0.5 


1.0 


7.81 


150 


Case 5: (Hyp.) locked Earth 


1.0 


1.0 


0.15 


5.23 


73 


Case 6: Small super-Earth (6 Mg) 


0.2 


0.5 


0.65 


4.15 


150 


Case 7: Large super-Earth (10 Me) 


0.2 


0.5 


0.96 


4.13 


150 


Case 8: Ocean Planet 


0.2 


0.5 


0.37 


2.80 


150 



Table 1 

Parameters for the different planetary configurations, d: orbital distance, M^: stellar 
mass, M.: planetary magnetic moment, R^. standoff distance, B: magnetic field at 
the magnetopause. 



distances (Section [2.3p . BotllfJ contribute to the magnetospheric compression, 
which in turn determines the flux of high energy cosmic ray particles onto the 
planetary atmosphere. Here, we compare the relative importance of these two 
effects. 



We first estimate the importance of the orbital distance via the stellar wind 
parameters. In Fig. [6l we first compare the cosmic ray energy spectra of the two 
strongly magnetised cases (with M. = I.OA^e): The case of the Earth (case 2, 
solid line) and the hypothetical case of a strongly magnetised exoplanet (case 

4, triangles). It can be seen that the cosmic ray energy spectra are identical 
at d =1.0 AU (case 2, solid line) and at d =0.2 AU (case 4, triangles). We 
also compare the two weakly magnetised cases (with A4 = O.ISA^e): the 
weakly magnetised exoplanet case with d =0.2 AU (case 3, dashed line) and 
the hypothetical case of the weakly magnetised Earth at d =1.0 AU (case 

5, circles). Again, despite the large difference in orbital distance (0.2 and 1.0 
AU) , Fig. E] shows that the cosmic ray energy spectrum is identical for these 
two cases. 



^ In principle, a third effect plays a role, namely the stellar mass and radius, which 
determine the stellar wind (e.g. via Eq. ([3])). However, if a solar mass star was used 
in case 3, the planetary magnetosphere would be even more compressed, which does 
not decrease the cosmic ray flux. Thus, the two other effects have to be responsible 
for the large difference in cosmic ray fluxes. 
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The distance-dependent stellar wind ram pressure does not seem to have any 
measurable influence on the cosmic ray energy spectrum. This may seem sur- 
prising, as it clearly has some influence on the size of the planetary magne- 
tosphere, see Tabled] (column 5). However, Tabled], column 6 shows that the 
magnetic fleld at the magnetopause is increased when the orbital distance of a 
planet is decreased while keeping the magnetic moment constant. This directly 
results from the pressure balance in eq. (jl]). Thus, when the orbital distance 
is decreased while the magnetic moment is kept flxed, the smaller size of the 
magnetosphere is compensated by a larger magnetic fleld at the magnetopause 
B{Rs), which keeps the cosmic ray shielding at a constant level. In other words, 
Rs, the size of the planetary magnetosphere, is not a sufficient indicator for 
the quality of the magnetospheric shielding of a planetary atmosphere. 

10 E — I 1 1 1 1 1 3 




proton energy [MeV] 

Fig. 6. Cosmic ray energy spectrum. Dash-dotted line: energy spectrum outside the 
magnetosphere (case 1, "unmagnetised case"). Dashed line: "locked exoplanet case" 
(case 3). Circles: hypothetical "locked Earth case" (case 5). Solid line: "(unlocked) 
Earth case" (case 2). Triangles: hypothetical "unlocked exoplanet case" (case 4). 

4-2 Influence of tidal locking 

If the difference of the cosmic ray energy spectra between a locked exoplanet 
and the Earth cannot be explained by the different orbital distances, the dis- 
crepancy lies in magnetic moment differences. This can be seen by comparing 
the cosmic ray energy spectrum of the Earth case (case 2, solid line) to the 
cosmic ray flux of a hypothetical Earth twin planet with a reduced magnetic 
moment (case 5, circles), with all other parameters kept constant {d = 1.0 AU, 
= I.OM0). Fig. [6] shows that the energy spectrum is much higher in the 
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weakly magnetised case {A4 = O.ISA^e, case 5, circles) than in the strongly 
magnetised case {A4 = I.OA^e, case 2, solid hne). We also compare the two 
exoplanet cases {d = 0.2 AU, = 0.5Mq), namely the "locked exoplanet 
case" (case 3, dashed line) and the "unlocked exoplanet case" (case 4, trian- 
gles). Again, it can be seen that the cosmic ray energy spectrum has much 
higher values for the weakly magnetised case than in the strongly magnetised 
case. Obviously, the influence of tidal locking, hence the reduced magnetic 
moment (and not the stellar wind ram pressure) is the decisive factor for the 
increased influx of cosmic ray particles in case 3 (when compared to case 2). 
The similarity of cases 2 and 4 (and of cases 3 and 5) also shows that protons 
of the chosen energy range mostly feel the planetary dip ole fleld, and that the 



fleld of the magnetpause currents has little influence (jVogt et al.l . 120071 . eq. 
20). 

Similarly to a change in orbital distance, the planetary magnetic moment 
directly influences the size of the planetary magnetosphere (see Table [1], col- 
umn 5). When the magnetic moment is changed, however, the magnetic fleld 
strength at the magnetopause (shown in Table [H column 6) does not dif- 
fer (compare cases 2 and 5, and cases 3 and 4, respectively). Thus, when 
B{Rs) =const, the quality of the magnetospheric shielding of the planetary 
atmospheres depends on the size of the magnetosphere Rs, i.e. it is weaker for 
those planets with a smaller magnetosphere. 



We conclude that, of the two effects suggested by iGrieBmeier et al.l (j2005bl ). 
only one (namely the effect of tidal locking) has a measurable influence on the 
cosmic ray flux. The other effect (stellar wind density variation with orbital 
distance) appears to have negligible influence. 



4-3 Influence of planetary type 



For planets with different size or composition than the Earth, the planetary 
magnetic dipole moment and the size of the magnetosphere are different from 
the Earth-like case. A few such cases are shown in Table [1] and discussed in 
this section. 

Both "super-Earth" cases (cases 6 and 7) correspond to the case of a terres- 
trial planet with a composition similar to that of the Earth, but with higher 
planetary masses (6 and 10 Earth masses, respectively) ar id correspondingly 
larger radii (1 . 63 an d 1.86 Earth radii, see the models of Leger et al.l . 20041 : 



Valencia et al.l . 120061 ). For both cases, the expected magnetic moment (Table 
[H column 4) is larger than for an Earth-like planet (case 3). The larger mag- 
netic moments lead to more extended magnetospheres (Tabled], column 5). 
Because the planetary radius is larger, however, the planetary atmosphere is 
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located at a larger distance from the planetary centre, so that the protection 
from galactic cosmic rays should not be expected to be much more efficient 
than for an Earth-like planet. Figure [7] shows that this is indeed the case: 
The energy spectra of the "locked exoplanet case" (case 3, dashed fine), of the 
"small super-Earth case" (case 6, double-dashed line) and of the "large super- 
Earth case" (case 7, triple-dashed line) are relatively similar. It can be seen 
that massive terrestrial planets do not necessarily have a much more efficient 
magnetospheric protection against galactic cosmic rays than an Earth-mass 
planet. 



The "Ocean Planet case" (case 8, dotted line) corresponds to a planet with a 
different composition: 17% metals, 33% sili cates, and 50 % wa ter ice (this cor- 
responds to the Ocean Planet modelled by lLeger et al.l . 12004 ). This change in 
composition reduces the average density when compared to the "super-Earth" 
case of the same mass: The radius of a 6 Earth-mass Ocean Planet is 2.0 Earth 
radii (as compared to the 1.63 Earth radii used above). This moves the atmo- 
sphere further away from the centre of the protecting dipole. Also, the reduced 
density leads to a somewhat reduced magnetic moment estimation. This is re- 
flected in the small standoff distance measured in planetary radii (see Table 
[T]) when compared to the "small super-Earth case" (case 6, double-dashed 
line), and the cosmic ray flux through the magnetosphere can be expected to 
be larger. These expectation is confirmed by the cosmic ray energy spectrum 
in Figure [71 More cosmic ray particles reach the atmosphere for the "Ocean 
Planet case" (case 8, dotted line) than that of the "small super-Earth case" 
(case 6, double-dashed lines). Note that both planets have the same mass (6 
Me), so that the difference is caused only by the difference in composition. For 
energies above 100 MeV, the cosmic ray protection of an Ocean Planet (case 
8) is weaker by a factor 2-3 than for a rocky planet of similar mass (case 6). 
We conclude that the planetary composition can be a potentially important 
factor concerning cosmic ray protection. 



5 Implications 



The increased flux of cosmic ray protons of galactic origin to the upper atmo- 
spheres of M-star planets has implications for the atmospheric composition 
and thus for the remote detection of biomarkers. It will also influence the 
flux of (secondary) cosmic ray particles reaching the ground, which can be 
expected to have biological implications. 
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proton energy [MeV] 



Fig. 7. Cosmic ray energy spectrum for different types of planets. Dasli-dotted line: 
energy spectrum outside the magnetosphere (case 1, "unmagnetised case"). Dotted 
line: "Ocean Planet case" (case 8). Dashed line: Earth-like planet (case 3, "locked 
exoplanet case"). Double-dashed line: "small super-Earth case" (case 6). Triple- 
dashed line: "large super-Earth case" (case 7). Solid line: "(unlocked) Earth case" 
(case 2). 

5.1 Atmospheric implications 



In the near future, atmospheres of extrasolar planets will be studied from the 
point of view of possible biological activity. "Atmospheric biomarkers" are 
compounds present in an atmosphere which imply the presence of life and 
which cannot be explained by inorganic chemistry alone. The simultaneous 
presence of significant amounts of atmospheric reducing gases (e.g. methane) 
and oxidising g ases (e.g. oxygen) as on the Earth is also an indicator of bio- 
logical activity (Lovelock, 1965 : Sagan et al. . 19931 ). Good biomarkers include 
oxygen (produced by photosynthesis), ozone (mainly produced from oxygen) 
and nitrous oxide (produced almost exclusively from bacteria). 



For the study of biomarkers in planetary atmospheres, it is important to know 
all inorganic effects wh ich can modify the abundances of these molec ules. 
Modelling studies (e.g. Schindler and Kasting . 2000l : Selsis et all . 2002 ) are 
required to rule out false c onclusions in cases whe re inorganic chemistry can 
mimic the presence of life. iGrenfell et al.l ( l2007bl ) looked at the response of 
biomarker chemistry on varying planetary parameters (e.g. orbital distance, 
stellar type). Here, we will discuss the response of biomarkers to cosmic rays. 



The influence of a high flux of OCRs to the atmospheres of planets in the 
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habitable zone (HZ) around M-stars was studied by iGrenfell et al.l (j2007al ). 
When cosmic rays travel through an Earth-like atmosphere, they have suffi- 
cient energy to break the strong N2 molecule and the chemical products react 
with oxygen to form nitrogen oxides (NOx). If UV levels are intense, e.g. cor- 
responding to the upper mesosphere and above on the Earth, then the NOx 
cannot survive. NOx affects biomarkers differently depending on altitude. In a 
tropo spheric environment , ozone is enhanced by NOx due to the smog mecha- 



nism (iHaagen-Smit 



19521). I n a st ratospheric environment, ozone is depleted 



by catalytic cycles (ICrutzenl . Il970l ). Once ozone is affected, the other biomark- 



ers can change too, because ozone affects the temperature profile (hence e.g. 
mixing processes, chemical sources and sinks) and UV levels (hence photolysis 
rates) of other chemical species. 



Using the cosmic ray spectra of our cases 1, 2 and 3, it was shown that GCRs 
may enhance the production of NOx by a factor of three, which may modify 
the a bundances of b i omarke r molecules (especially water and ozone) by up to 
25% fiGrenfell et al.l . l2007ah . 



It thus seems that most biomarkers are not strongly influenced by galactic 
cosmic rays for Earth-like planets in close orbits around M-stars. The case is 
different for solar energetic particles, w hich can strongly mod ify the abundance 
of ozone in the planetary atmosphere (IGrenfell et al.l . l2008l ). Such changes in 
atmospheric biomarker concentrations have to be taken into account when 
searc hing for biosignatures in the spectra of Earth-like ex oplanets (jSelsid . 
2004h bv future space rn issions like DARWIN dPridlundl . liool or SEE-COAST 
( ISchneider et al.l . 120061 ). 



5.2 Biological implications 



While a quantitative treatment is not yet available, the increased flux of cosmic 
rays at the top of the planetary atmosphere has not only implications for atmo- 
spheric chemistry and biomarkers, but also for the flux of secondary radiation 
on the surface. When galactic cosmic rays of sufficiently high energy reach 
the planetary atmosphere, they generate showers of secondary cosmic rays. Of 
the different components of such showers, slow neutrons have the strongest 
influence on biological systems. The propagation of high energy cosmic ray 
particles also depends on the composition and density of the planetary atmo- 
sphere. For an Earth-like atmosphere, the minimum energy which a proton 
must have to initiate a nucl ear interaction that may be detectable at se a level 
is approximately 450 MeV ( Reeves et al. . 1992 : Shea and Smart , 2000l ). Ter- 
restrial planets with dense atmospheres like Venus (100 bar surface pressure) 
would be better shielded by the planetary atmosphere, so that no secondary 
radiation can reach the surface. On the other hand, for planets with thin at- 
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mospheres like Mars (6 mbar surface pressure), the surface would probably be 
totally sterilised, even for a relatively low cosmic ray flux. 



Keeping the atmospheric pressure constant at ~1 bar, the flux at energies 
> 450 MeV is strongly increased for the planetary situations presented in 
Section m (see Figs. [6] and [7]), and a large increase of secondary cosmic rays 
can be expected at the planetary surface. 

Biological effects, namely an increase in cell fusion indices for different cell- 
hnes, were found to be significantly correlated with the neut r on cou nt rate at 



the Earth's surface (IBelisheva et al.l . l2005l : iGrieBmeier et al.l . l2005bl ). Similar 



but much stronger and more diverse effects were observed during large solar 
particle events, where solar cosmic rays dominate over galactic cosmic rays. 
The effect of cosmic rays on cells and on DNA are potentially hazardous for 
life. On the other hand, changes at the genetic level are a necessary condition 
for biological evolution, so that cosmic rays may also play a role in radiation- 
induced evolutionary events. 

As tidally locked Earth-like exoplanets inside the HZ of low-mass K / high- 
mass M stars are only weakly protected against high energetic cosmic rays, 
they can be expected to experience a higher surface neutron flux and stronger 
biological effects than Earth-like planets with a strong magnetic moment. For 
this reason, it may be more difficult for life to develop on the surface of planets 
around such stars. It was also shown (Section 14. 3p that larger planets (super- 
Earths) are not much better protected against GCRs than smaller planets. 



6 Conclusion 



The atmospheres of extrasolar planets with orbits in the habitable zone of M 
stars are exposed to high fluxes of galactic cosmic rays. We have studied the ef- 
fect of orbital distance via the stellar wind parameters, tidal locking, planetary 
mass and planetary composition on the cosmic ray fluxes. In previous works, 
the weak shielding found for M star planets was attributed to a combination of 
the enhanced ste llar wind ram pre s sure w ith the reduced planetary magnetic 
dipole moment ( Griefimeier et al. . 2005bl lal). Here, we have shown that the 



variable stellar wind density at different orbital distance ranging from 0.2 to 
1.0 AU (i.e. the GHZ of star with O.SM© < < l.OM©). has little influence 
on the GGR flux. The planetary magnetic moment is the decisive parame- 
ter: Tidally locked planets with a small magnetic moment have much higher 
particle fluxes than freely rotating planets with stronger magnetic moments. 

We also find that more massive planets are not necessarily much better pro- 
tected against galactic cosmic rays, but that the planetary composition can 
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play an important role. 



Generally speaking, to determine the magnetosplieric shielding of a planetary 
atmosphere, it is not sufficient to calculate Rg, the size of the planetary mag- 
netosphere, but one has also to take into account B{Rs), the magnetic field at 
the magnetopause. 

The above factors have to be considered for missions studying biosignatures 
in the spectra of Earth-like exoplanets for two reasons: The enhanced flux 
of galactic cosmic rays m ay naturally generate or destroy molecules typically 



regarded as biomarkers (iGrenfell et al.l . l2007al ). which has to be taken into 



account in the analysis and interpretation of such spectra. Also, as the habit- 
ability of M star planets may be affected, future observation campaigns should 
ensure that enough non-M stars are in the sample to maximise the likelihood 
of a positive detection. 
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